DEPARTMENT  OF  THE  NAVY 
NAVAL  SHIP  RESEARCH  AND  DEVELOPMENT  CENTER 

Betteida,  Md.  200V. 


5 


SIMPLIFIED  THEORETICAL  METHODS  OF  PREDICTING 
THE  MOTIONS  OF  A  CATAMARAN  IN  WAVES 


by 

R.  Wahab  and 
E.N.  Hubble 


I 


APPROVED  FOR  PUBLIC  RELEASE:  CISTRIBUIION  UNLIMITED 


I 


January  1972 


Report  3736 


TABLE  OF  CONTENTS 


Page 

ABSTRACT  .  1 

ADMINISTRATIVE  INFORMATION  . . .  1 

INTRODUCTION  . .  1 

MOTION  PREDICTION  METHODS  .  2 

PITCH  AND  HEAVE  . 2 

ROLL  .  4 

COMPARISON  OF  COMPUTED  AND  MEASURED  DATA  .  7 

HEAVE  . 7 

Theory  versus  Experiment  for  the  Catamaran  . . .  7 

Catamaran  versus  Conventional  Ship  . 11 

PITCH  .  11 

Theory  versus  Experiment  for  the  Catamaran  .  11 

Catamaran  versus  Conventional  Ship  . . .  11 

ROLL  .  11 

CONCLUDING  REMARKS  .  12 

APPENDIX  A  -  PROCEDURE  AND  NOTATION  USED  IN  COMPUTER  PROGRAM  RLAC 
FOR  PREDICTING  THE  ROLL  OF  A  CATAMARAN  IN  REGULAR  AND 
IRREGULAR  SEAS  .  13 

APPENDIX  B  -  FORMAT  OF  INPUT  FOR  PROGRAM  RLAC  . .  19 

APPENDIX  C  -  SAMPLF  OUTPUT  FROM  PROGRAM  RLAC  . 23 

APPENDIX  D  -  FORTRAN  LISTING  OF  PROGRAM  RLAC  . 29 

REFERENCES  . 35 

LIST  OF  FIGURES 

Page 

Figure  1  -  Hull  Lines  of  Catamaran  Model  5061  .  3 

Figure  2  -  Lewis  Sections  with  the  Same  Waterline  Width,  Draft,  and 

Sectional  Area  as  One  Hull  of  the  Catamaran  .  3 

Figure  3  -  Lewis  Sections  with  the  Same  Waterline  Width,  Draft,  and 

Sectional  Area  as  Both  Catamaran  Hulls  .  5 

Figure  4  -  Particulars  of  Catamaran  Model  3061  . .  ,  8 

Figure  5  -  Metacentric  Height  and  Transverse  Gyradius  of 

Model  5061  8 

Figure  6  -  Heave  Transfer  for  Head  Waves  and  Various  Speeds  .  ? 

Figure  7  -  Pitch  Transfer  for  Head  Waves  and  Various  Speeds  .  . .  9 

Figure  8  -  Roll  Transfer  in  Beam  Waves  for  Various  Hull 

Separations  . . .  •  10 

ii 


NOTATION 


Area  of  waterplane  of  one  hull 
B  Breadth  of  one  hull 

B^  Transverse  distance  from  longitudinal  axis  of  symmetry  of  catamaran 
to  longitudinal  axis  of  one  hull 

B  Overall  breadth  of  catamaran 

m 

BG  Vertical  distance  of  center  of  gravity  above  center  of  buoyancy 
b  Coefficient  of  the  restoring  moment 

CG  Center  of  gravity 

F  Exciting  moment 

g  Gravitational  acceleration 

GM  Transverse  metacentric  height 

I  Moment  of  inertia  of  waterplane  area  of  both  hulls  with  respect  to 

°  to  the  longitudinal  axis  of  symmetry  of  the  catamaran 

I_  Moment  of  inertia  of  waterplane  area  of  one  hull  with  respect  to  its 

longitudinal  axis  ' 

Transverse  gyradius 

L  Hull  length  between  perpendiculars 

m  Mass  moment  of  inertia 

n  Damping  coefficient 

t  Time 

V  Speed  of  advance 

zA  Heave  amplitude 

e  Phase  angle 

5a  Wave  amplitude 

^  Wave  height 

0A  Pitch  amplitude 

tc  Wave  number  =  2tt/X 

X  Wavelength 

p  Water  density 

d>A  Roll  amplitude 

u  Circular  frequency 


in 


Vj  Volume  of  water  displaced  by  one  hu!l 

Volume  of  water  displaced  by  both  catamaran  hulls 


ABSTRACT 


Simplified  methods  are  discussed  for  estimating  (1)  the 
pitch  and  heave  of  catamarans  in  head  seas  based  on  theory  which 
has  proven  successful  for  conventional  single  hulled  ships,  and 
(2)  the  roll  of  catamarans  in  beam  seas  by  representing  the 
small  amount  of  roll  as  alternate  heaving  of  the  two  hulls. 

Both  prediction  methods  neglect  interaction  effects  between 
the  two  hulls.  Computed  values  of  pitch,  heave,  and  roll  are 
compared  with  experimental  data  from  model  tests  of  a  catamaran 
in  regular  waves.  Documentation  of  the  computer  program  for 
predicting  the  roll  of  a  catamaran  in  regular  and  irregular  seas 
is  presented  in  the  appendices . 

ADMINISTRATIVE  INFORMATION 


This  work  was  performed  at  Naval  Ship  Research  and  Development 
Center  (N3RDC)  primarily  under  the  Naval  Ship  Systems  Command  (NAVSHIPS) 
Exploratory  Development  Applied  Hydromechanics  Program,  Subproject  SF  35. 
421.006,  Task  1713.  Development  of  the  computer  routine  for  predicting 
roll  in  beam  seas  was  undertaken  as  part  of  a  conceptual  research  feasibi¬ 
lity  study  of  catamaran  aircraft  carriers  and  funded  from  NSRDC  in-house 
Project  1-H71-001 ,  Task  2F  35.412.002. 

INTRODUCTION 

The  growing  interest  in  catamarans  makes  it  desirable  to  be  able  to 
predict  the  motions  of  these  ships  by  techniques  similar  to  those  which 
have  been  developed  for  monohulls.  Existing  computer  programs  for  pre¬ 
dicting  the  pitch  and  heave  motions  of  single-hulled  ships  provide  a  first 
approach  for  predicting  the  pitch  and  heave  of  catamarans  in  head  seas. 

The  basic  assumption  in  the  present  approach  is  that  the  hulls  are  widely 
separated,  i.e.,  interaction  effects  between  the  two  hulls  are  neglected. 
With  this  assumption,  it  is  relatively  simple  to  write  a  computer  program 
for  estimating  the  roll  motion  of  a  catamaran  in  beam  seas.  Since  rolling 
of  a  catamaran  takes  place  with  small  angles,  it  can  be  regarded  as  alter¬ 
nate  heaving  of  the  two  hulls,  therefore,  parts  of  the  program  to  compute 
pitch  and  heave  can  be  used  for  the  prediction  of  roll  of  a  catamaran.  In 
thi •  report  the  motions  estimated  in  the  manner  described  above  are  com¬ 
pared  with  experimentally  obtained  data  from  catamaran  Model  5061  which  has 
been  tested  at  this  Center  with  various  hull  separations. 
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Since  the  hulls  of  a  catamaran  generally  have  proportions  different 
from  those  of  a  conventional  ship  hull,  the  effect  of  the  beam-draft  ratio 
on  the  motions  in  head  waves  is  also  examined  to  some  extent. 

MOTION  PREDICTION  METHODS 


PITCH  AND  HEAVE 

The  pitch  and  heave  motions  of  conventional  ships  in  head  waves  at 
Froude  numbers  up  to  0.45  have  been  predicted  quite  successfully  using  the 
Frank  Close-Fit  Ship-Motion  Computer  Program  YF17.*  The  regular  wave 
responses  are  computed  according  to  an  improved  version  of  the  Korvin- 
Kroukovsky  strip  theory.  An  essential  part  of  the  program  is  the  compu¬ 
tation  of  the  sectional  added  mass  and  damping  coefficients  by  either  the 
Lewis-form  method  or  the  more  accurate  but  time-consuming  close-fit  method. 
The  same  program  (hereafter  referred  to  as  YF17)  has  been  used  for  the 
calculation  cf  catamaran  pitch  and  heave  in  head  seas  presented  in  this 
report.  The  catamaran  considered  here.  Model  5061,  has  hulls  with 
asymmetric  sections  forward  of  midship;  see  Figure  1.  However,  YF17  con¬ 
siders  only  a  single  body  which  is  symmetrical  about  a  vertical  longi¬ 
tudinal  plane.  Therefore,  in  the  equations  of  motion  the  added  mass  and 
the  damping  coefficient  computed  for  each  section  were  those  for  a  Lewis 
section  having  the  same  waterline  width,  draft,  and  sectional  area  as  one 
hull  of  the  catamaran.  These  sections  are  shown  in  Figure  2.  This  Lewis- 
form  method  has  been  used  in  lieu  oZ  the  close-fit  method  for  many  con¬ 
ventional  ships  (except  those  with  large  bulbous  bows)  without  significantly 
effecting  the  resultant  computed  motions.  Experience  gained  with  compu¬ 
tation  procedures  which  differ  only  slightly  from  those  used  in  YF17,  in 
combination  with  Lewis  sections,  indicates  that  the  agreement  between  ex¬ 
periment  and  theory  is  better  for  beamy  hulls  than  for  rather  nairow  hulls; 

2  3 

see  Jooseu  et  al.  and  Vassilopoulos  and  Mandel.  Catamaran  hulls 
generally  belong  to  the  latter  category;  those  considered  in  this  report 


References  are  listed  on  page  33. 
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Hull  Lines  of  Catamaran  Model  5061 
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have  a  beam-draft  ratio  of  1.3.  To  investigate  the  effect  of  the  beam- 
draft  ratio,  the  computed  motions  in  head  waves  of  the  catamaran  with  the 
narrow  hulls  are  compared  with  the  computed  motions  of  a  ship  that  has  the 
same  length,  draft,  and  displacement  as  the  catamaran,  and  waterline 
widths  and  sectional  areas  equal  to  those  of  both  catamaran  hulls.  The 
sections  of  the  conventional  ship  are  given  in  Figure  3.  Because  of  the 
limitations  of  the  prediction  method  for  catamarans,  the  compared  motions 
of  the  two  ships  cannot  be  regarded  as  correct  in  the  quantitative  sense, 
but  only  qualitatively. 

ROLL 

A  slight  modification  of  the  theory  outlined  by  Wahab^  was  used  to 
develop  a  computer  program  for  predicting  the  rolling  characteristics  of  a 
catamaran  in  both  regular  and  irregular  seas.  Complete  documentation  for 
this  program,  designated  RLAC,  is  presented  in  Appendixes  A-D. 

The  theory  is  based  on  the  assumption  that  the  rolling  of  a  cata¬ 
maran  can  be  represented  by  alternate  heaving  of  the  two  hulls  without 
significant  error  since  the  roll  angles  as  well  as  the  roll  damping  and 
added  moments  of  inertia  of  each  hull  are  small. 

In  determining  the  exciting  moment  in  beam  waves,  it  was  assumed 
that  the  presence  of  the  ship  did  not  change  the  pressure  distribution  in 
the  undisturbed  wave.  The  exciting  moment  was  obtained  from  the  hydro¬ 
static  pressure  acting  on  the  ship  with  a  correction  for  the  Smith  effect. 
This  approach  is  known  to  give  reasonable  results  in  head  waves ,  but  no 
verification  has  been  made  for  the  case  of  beam  waves. 

The  uncoupled  linear  equation  of  motion  is 

•  _ 

m  (j)  +  n<J>  +  b<J>  =  F  sin  cot  (1) 

After  the  starting  transient  has  died  out,  the  solution  of  this  equation  is 

=  <p^  sin  (ait  +  e)  (2) 
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=  Fy/j/cb  -  m  w2)2  +  n2  u ? 


(3) 


u>  n 


e  =  atan 


m  ij)  -  b 


(4) 


For  regular  beam  waves  the  exciting  moment  F  is 


F  -  BcAWe<' 


17V  sin  (<  Bc)*  2  +  k  (I  -  BG  cos  (k  Bc)  ‘  (5) 


where  k  -  2ir/A. 


The  coefficient  b  of  the  restoring  moment  may  be  ~alculated  by 


b  =  GM  pg  V2  =  (Io/72  -  BG)  pg  V2 


(6) 


When  the  catamaran  rolls  with  amplitude  <j)^,  each  hull  heaves  with  amplitude 
<f>ABc  in  addition  to  the  rolling.  Therefor^,  the  mass  moment  of  inertia  is 
subdivided  as  follows: 


m  =  m  +  m, ,  +  B  “m 
c  cpd>  c  zz 


(7) 


where  m  is  the  transverse  moment  of  inertia  of  the  catamaran  itself, 
c 

m^  is  the  added  moment  of  inertia  due  to  rolling  of  both  hulls,  and 
mzz  is  the  added  mass  due  to  the  heaving  motion  of  both  hulls. 

2 

Since  m^  is  small  compared  to  Bc  mzz>  if  is  neglected. 

The  damping  coefficient  can  be  subdivided  as  follows: 


n  -  B  n  +  n, . 
c  zz  ®<p 


(8) 


where  n^  is  the  damping  coefficient  due  to  rolling  motion  of  both  hulls 

and  n  is  the  damping  coefficient  due  to  heaving  motion  of  both  hulls, 
zz  2 

Since  n^  is  small  compared  to  B£  nzz,  it  is  also  neglected. 

Program  RLAC  incorporates  Subroutines  ADMAB  and  NILS  from  Program 
YF17  for  computing  the  added  mass  and  damping  coefficient  due  to  heave.  As 
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with  pitch  and  heave,  the  Lewis -form  sections  shown  in  Figure  2  are  used 
for  the  roll  computation  of  the  catamaran  being  studied,  and  interaction 
effects  between  the  two  hulls  are  neglected. 

In  view  of  the  aforementioned  limitations  of  the  existing  theory  of 
catamaran  roll,  refinements  such  as  (1)  correction  for  forward  speed 
effects  on  the  coefficients  of  the  equations  of  motion  and  (2)  correction 
to  the  exciting  moment  for  added  mass  and  damping  forces  associated  with 
the  oribital  motion  of  the  water  particles  in  the  waves  have  not  been 
made. 

Program  RLAC  can  also  be  used  to  predict  catamaran  roll  in  irregular 
seas.  Roll  is  computed  for  a  range  of  wave  frequencies  using  an  arbitrary 
wave  steepness  (?WA)  of  1/50.  The  method  of  linear  superposition  on  the 
sea  spectrum  given  by  the  Pierson-Moskowitz  formulation  is  used  for  pre¬ 
diction  of  the  roll  displacement  and  acceleration  at  various  sea  states. 
Calculations  for  catamaran  Model  506’  at  significant  wave  heights  of  4,  10, 
20,  and  30  ft  are  contained  in  the  sample  output  shown  in  Appendix  C. 

COMPARISON  OF  COMPUTED  AND  MEASURED  DATA 

The  particulars  of  catamaran  Model  5061  and  the  dynamic  conditions 
for  which  it  was  tested  are  given  in  Figures  4  and  5,  respectively.  In 
Figures  6  and  7  the  computed  heave  and  pitch  motions  are  compared  with 
results  of  experiments  in  head  waves.  The  dashed  curves  in  the  figures 
represent  the  computed  motions  of  a  ship  with  the  same  length,  draft,  and 
displacement  as  the  catamaran  and  with  waterline  widths  and  sectional  areas 
equal  to  those  of  both  of  the  catamaran  hulls.  A  comparison  between 
computed  and  measured  roll  motions  for  the  catamaran  is  made  in  Figure  3. 

HEAVE 

Theory  versus  Experiment  for  the  Catamaran 

For  zero  speed  and  all  hull  separations,  it  is  seen  in  Figure  6 
that  the  computed  values  agreed  well  with  measured  Leave  except  for 
A/L  =  1.1  where  a  slight  peak  was  obtained.  Trends  were  maintained  for  the 
remaining  speeds.  Ho'ever,  measured  amplitudes,  especially  in  vhc 
resonance  regior,  were  significantly  lower  than  predicted  for  A  >  L  but 
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Figure  4  -  Particulars  of  Catamaran  Model  5061 


Transfer  for  Head  Waves  and  Figure  7  -  Pitch  Transfer  for  Head  Waves 

Various  Speeds  Various  Speeds 


•«%# 


was  slightly  larger  than  predicted  for  X  <  L  in  most  cases.  It  has  been 
2  3 

shown  *  that  the  present  state-of-the-art  calculation  procedure  yields 

poor  results  even  for  conventional  ships  with  low  beam-draft  ratios.  It 

overestimates  the  pitch  ar.d  heave  response  amplitudes,  particularly  at 

3 

resonance.  Vassilopoulos  and  Mandel  attribute  this  to  the  use  of  Lewis 
sections.  The  discrepancies  found  in  the  present  comparison  may  possibly 
also  be  partly  attributed  to  imperfections  in  the  theory  as  applied  to 
catamarans . 

Catamaran  versus  Conventional  Ship 

The  curves  in  Figure  6  indicate  that  a  catamaran  may  be  expected  to 
heave  more  than  a  monohull  ship  with  the  same  length,  draft,  and  displace¬ 
ment.  However,  since  the  theory  overestimates  the  motions  for  low  beam- 
draft  ratio  hulls,  the  difference  will  actually  be  smaller  than  the  two 
computed  curves  indicate. 

PITCH 

Theory  versus  Experiments  for  the  Catamaran 

For  the  two  lowest  investigated  speeds,  the  experimentally  obtained 
pitch  shown  in  Figure  7  had  about  the  same  trend  as  the  predictions.  For 
the  two  higher  speeds,  however,  there  was  a  distinct  difference  in  the 
character  of  the  measured  and  computed  transfer  curves;  this  indicates  that 
the  interaction  effects  between  the  hulls  are  most  likely  not  negligible. 

To  some  extent,  the  difference  may  also  be  due  to  imperfections  in  the 
computation  procedure  as  discussed  in  the  previous  section. 

Catamaran  versus  Conventional  Ship 

The  pitch  motion  of  a  catamaran  may  be  expected  to  be  large  com¬ 
pared  to  that  of  a  ship  with  the  same  length,  draft,  and  displacement,  be¬ 
cause  of  the  smaller  beam-draft  ratio  of  the  catamaran  hulls. 

ROLL 

The  general  nature  of  the  roll  behavior  in  Figure  8,  i.e.,  slope 
and  location  of  maxima,  agreed  fairly  well  with  prediction.  It  is  noted 
that  the  computation  is  primarily  valid  for  zero  speed  since  no 
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spe  ;d-dependent  terms  were  included  in  the  equation  of  motion.  However,  on 
the  basis  of  the  experimental  data,  some  functional  dependence  of  roll 
damp.-.g  on  hull  separation,  and  to  a  lesser  extent  on  forward  speed,  is 
apparent.  Therefore,  it  may  be  worthwhile  to  refine  the  equation  of  motion 
but  maintain  the  assumption  of  no  mutual  influence  be’iveen  the  hulls. 

CONCLUDING  REMARKS 

it  appears  that  neglecting  the  interaction  effects  between  the  hulls 
does  not  prevent  reasonable  results  when  computing  roll  in  beam  seas. 

Better  results  may  possibly  be  obtained  by  including  speed-dependent  terms 
in  the  ,  quation  of  motion.  The  pitch  and  heave  motions  in  head  waves  could 
not  be  satisfactorily  predicted  et  the  high  Froudc  numbers  of  0.25  and  0.38. 
The  discrepancy  may  be  partly  attributed  to  the  unsatisfactory  performance 
of  the  calculation  procedure  for  low  beam-draft  ratio  hulls. 

The  computations  also  showed  that  because  of  the  small  beam-draft 
ratio  of  its  hulls,  the  behavior  of  a  catamaran  in  head  waves  may  be  sig¬ 
nificantly  worse  than  the  behavior  of  a  ship  with  the  same  length,  draft, 
and  displacement. 
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APPENDIX  A 

PROCEDURE  AND  NOTATION  USED  IN  COMPUTER  PROGRAM  RLAC 


BASIC  HULL  GEOMETRY 


*  BPL  =  Lgp  =  length  between  perpendiculars  of  each  hull  (ft) 

*  NOS  =  n  =  number  of  stations 

s 

*  ST(K)  =  Sta.  =  station  number  (Sta  0  must  be  at  the  FP) 

(Sta  20  must  be  at  the  AP) 

*  NM(K)  =  =  number  of  waterlines  at  which  offsets  are  given 

X(K)  =  =  distance  of  Sta^  aft  of  FP  (ft)  =  Sta^‘Lgp/20 

lf  \  * 0 

•» 

*  B(K)  =  b^  =  full  beam  of  one  hull,  at  the  waterline  (ft) 

*  H(K)  =  h^  =  distance  from  keel  to  waterline  (ft) 

*  CA(K)  =  C.  =  area  coefficient  I 

AR(K)  =  A^  =  sectional  area  (ft  )  =  cA^*bjc'hjc  J 

If  >  0 

*  Z(J,K)  =  z.  ,  =  distance  above  the  baseline  (ft)  \  | 

-  ’  (Zj  ^  be  at  the  keel)  1  I 

(z  ,k  must  be  at  the  waterline  ^,mk^| 

“k  IS 


(k=l,ns) 


l,ns) 


Y(J,K) 

=  /j,k 

=  half  beam 

B(K) 

=  bk 

‘  2’>yk 

H(K) 

=  hk 

=  \-k  • 1 
* 

AR(K) 

=  Ak 

=2  J  y  dz 

CA(K) 

n 

=  V'V’v 

t 

Input  values. 

(k=l,ns) 


(numerical  intergration  by 
the  Simpson  rule) 


Note:  The  FORTRAN  designation  for  the  variables  is  given  in  the  first 
coltmm,  and  the  normal  notation  in  the  second  column. 
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AM(K)  =  A^ 

MS  =Q5 
EM  -  B 

HM  =  H 

RKO  ^  p 

G  =  g 

VOL1  =  V 

VOL  =  V2 

TM  =  M 

DLBS  =  A 

DTONS  = 

AW  =  Aw 

OIP  =  iT 

CB  =  C_ 

B 

CW  =  cw 

BOY  =  LCB 

CBL  =  LCB/LBp 

FLC  =  LCF 

CFL  =  LCF/Lgp 

BL  =  L/B 

BT  =  B/H 


moment  of  the  area  of  Staj,  about  the  FP  (k=l,ns) 
value  of  k  where  Sta^  =  1C 

full  bean  of  one  hull  at  amidships  (ft)  =  b^ 
draft  (keel  to  WL)  at  amidships  (>)  =  hjg 
water  density  =  1.9905  lb-sec^/ft^ 

2 

acceleration  of  gravity  =  32.174  ft/sec 

3  fL 

volume  of  water  displaced  by  one  hull  (ft  )  =  I 

3 

volume  displaced  by  both  hulls  (ft  )  =  2*V 

2 

total  mass  of  the  catamaran  (lb-sec  /ft)  =  pV^ 
displacement  of  catamaran  (lb)  =  p  g  V2 
displacement  of  catamaran  (tons)  =  A/2240 

2  fL 

area  of  waterplane  of  one  hull  (ft  )  =  I  b  dx 

J0 

moment  of  inertia  of  waterplane  area  of  one  hull 
with  respect  to  the  longitudinal  axis  of  the  hr  11 

=  2/3  f  b3  dx 
J0 

block  coefficient  of  one  hull  =  V^/(Lgp*B*H) 
waterplane  coefficient  of  one  hull  =  A^/(Lgp*B) 
distance  of  center  of  buoyancy  aft  of  FP  (ft) 


f  A  x  dx  /  f  A  dx 

-  0  J  L*'o  -I 


distance  of  center  of  floatation  aft  of  FP  (ft) 
LCB  +  f"  f  (x-LCB)  b  dxl  /  [A  1  =  f  f  b  x  dx 


r  r L  *1 

rL  i  / 

U  (x-LCB)  b  dx J  ^ 

/  [AJ  -  _ 

J0  bxdxJ/[A»1 

V8 


OTHER  SHIP  PARAMETERS 


BP 

L 

CL 

CG 

CB 

DK 

= 

KD 

GK 

= 

KG 

BK 

= 

KB 

GD 

= 

GD 

BG 

= 

BG 

GM 

= 

GM 

YL 

= 

B 

c 

YLP 

= 

Bd 

RG 

= 

k» 

01 

= 

I 

o 

CRM 

= 

b 

m 

c 

LING  MOTIONS  IN 

H21 

- 

V* 

WS 

= 

length  between  perpendiculars  of  each  hull  (ft) 

Note:  If  this  length  differs  from  the  one  used 
for  the  basic  hull  geometry  calculations,  all 
the  basic  parameters  (V,  B,  H,  etc.)  are  scaled 
by  the  appropriate  linear  ratio. 

centerline 

center  of  gravity 

center  of  buoyancy 

vertical  distance  of  deck  above  keel  (ft) 

vertical  distance  of  CG  above  keel  (ft) 

vertical  distance  of  CB  above  keel  (ft) 

distance  of  deck  above  CG  (ft)  =  KD  -  KG 

distance  of  CG  above  CB  (ft)  =  KG  -  KB 

metacentric  height  (ft)  =  KB  +  BM  -  KG  =  Ip/^  - 

transverse  distance  from  CL  of  catamaran  to 
CL  of  one  hull  (ft) 

transverse  distance  from  CL  of  catamaran  to 
outer  edge  of  the  deck  (ft) 

transverse  gyradius  (ft) 

I 

moment  of  inertia  of  the  waterplane  area  of  both 
hulls  with  respect  to  the  longitudinal  axis  cf 

2 

symmetry  of  the  catamaran  =  2  (1^,  +  B^  A^) 
coefficient  of  the  restoring  moment  =  GM  A 

transverse  moment  of  inertia  of  the  catamaran 
-  *  *  P  V, 


wave  height  to  length  ratio  =  1/50 
wave  slope  =  tt/50 


* 

Input  values. 


BG 
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*  OMIN 

= 

= 

minimum  nondimens ional  wave  frequency, 
generally  0.2 

*  OMAX 

— 

Wi)-, 

* 

maximum  nondimens ional  wave  frequency, 
generally  10.0 

*  DOM 

= 

■(4) 

= 

increment  of  nondimensional  frequency, 
generally  0.2 

NFR 

= 

nF 

= 

number  of  frequencies  =  [ (OMAX-OMIN) /DOM]  +  1 

Calculated 

for 

each  of  the 

nF 

frequencies : 

OMLG(N) 

- 

s 

*  4 Vi) 

OM(N) 

- 

U> 

= 

wave  frequency  (rad/ sec)  =  j 

WL 

= 

A 

s 

2 

wavelength  (ft)  =  2  tt  g/w 

WLL 

= 

A/L 

= 

ratio  of  wavelength  to  ship  length 

WH2 

s 

= 

wave  height  (ft)  =  X/50 

WH 

= 

CA 

= 

wave  amplitude  (ft) 

**  A33(N) 

- 

a33 

= 

added  mass  due  to  heave  of  each  hull/fpVp 

**  B33(N) 

= 

b33 

= 

damping  coefficient  for  each  hull  / (pV^  ^g/L) 

m 

zz 

= 

added  mass  due  to  heave  of  both  hulls 
=  a3J  PV2 

n 

zz 

= 

damping  coefficient  due  to  heave  of  both 
hulls  =  b33  pV2  y'g/L 

★ 

Input  values. 

**Values  of  a^7  and  b33  may  be  input  or  calculated  in  the  program.  If 

calculated  by  this  program,  the  sections  are  represented  by  the  Lewis-form 
method,  and  the  two-dimensional  added  mass  and  damping  coefficients  are 
calculated  according  to  the  Grim  method  by  Subroutine  ADMAB,  which  is  ab¬ 
stracted  from  Program  YF17,  but  was  initially  written  by  Stevens  Institute 
of  Techrology.  If  a33  and  b3-  are  input  directly,  they  may  be  obtained 

from  Program  YF17  which  uses  either  the  Lewis-form  or  the  close-fit  method 
for  each  section  independently,  as  desired.  In  either  case,  the  three- 
dimensional  values  are  computed  according  to  strip  theory  by  using  Sub¬ 
routine  NILS  (also  abstracted  from  YF17)  for  computation  of  the  Simpson 
weight  coefficients. 
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CM 

jr 

a 

a 

mass  moment  of  inertia  =  mc  ♦ 

CN 

s 

n 

3 

2 

damping  coefficient  =  B  n 

c  zz 

PWL 

3 

K 

= 

wave  number  =  2tt/a 

ANG 

= 

K  B 

c 

3 

2itB  /X 
c 

FBAR 

3 

3 

exciting  moment 

PHIB  = 


PHI  = 


RAOR  = 


RAOA  = 


<  C. 


(*)' 


3  2  m 
c  zz 


Cwpg 


^7 


/  ? 2 
ftr)  ' 


2  2 
K  V  /A  1  I  “ 

1  W  sin  (<Bc)  +  k(It-BG  Vp  cos  (tc  B£)| 


>] 


2,2  .  2  2 


roll  amplitude  (rad)  =  F  /  ]/{b  -  mu)t)''  +  n"  w 
roll  amplitude  /  wave  slope 


response  amplitude  operator  for  roll  displacement 
(rad/ ft) 2 


response  amplitude  operator  for  roll  acceleration 
(rad/ft/sec2)2 


ROLLING  MOTIONS  IN  IRREGULAR  WAVES 


*  NSWH 

nH 

=  number  of  significant  wave  heights  for 
irregular  sea  computations 

*  H13(M) 

Hl/3 

m 

=  significant  wave  height  j 

=  average  of  the  highest  one-third  wave 
heights  | 

(m=l,nH) 

Calculated 

for  each 

Hl/3’ 

to  combination: 

SW(N) 

s(o>) 

=  Pierson-Moskowitz  sea  spectral  formulation 
(ft2  sec) 

=  (A/a )^)  ,  where  A  =  0.0081  g2  and 

B  =  33.56/ (H1/3) 2 


Input  values. 
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FROLL (N,M) 
FACC (N,M) 


(4>a/^a)2  s(ui)  (rad2  sec) 

(V2/^a)2  s(u)  (rad2/sec3) 

Calculated  for  each 

Ei  = 2  /  (V’A)2  s(u)  **  (rad2) 

=  amplitude  of  significant  roll  angle  (rad)  = 

1.41 

t 

RDEG(M)  =  4»1/_  =  amplitude  of  significant  roll  angle  (deg)  = 

'  4>1/3  *  (180/it) 

E2  =  2  J  (4>aw2/£a)2  s(w)  d^  (rad2/ sec4) 

a^yj  =  amplitude  of  significant  roll  acceleration 
(rad/sec2)  =  1.41 

V 

ADEG(M)  =  ai/-$  =  amplitude  of  significant  roll  acceleration 

(deg/sec2)  =  a^3  (180/rr) 

AVG(M)  =  a  =  (amplitude  of  significant  vertical  roll 

vl/3  acceleration  at  outer  edge  of  deck)/ (gravi¬ 
tational  acceleration)  =  a^3  B^/g 

AHG(M)  =  a^  =  (amplitude  of  significant  horizontal  roll 

nl/3  acceleration  on  the  deck)/ (gravitational 

acceleration)  =  a ^  GD/g 
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APPENDIX  B 

FORMAT  OF  INPUT  FOR  PROGRAM  RLAC 


CARD  SET  1  (one  card) 


COLUMNS 

FORMAT 

FORTRAN 

2-72 

12A6 

TITLE 

CARD  SET  2  (one  card) 

COLUMNS 

FORMAT 

FORTRAN 

1-  9 

F9.3 

BPL 

10-18 

F9.3 

0MIN 

19-27 

F9.3 

0MAX 

28-36 

F9.3 

DOM 

37-45 

F9.3 

CST 

46-54 

19 

NOS 

55  -63 

19 

IAMD 

CARD  SET  3  (one  card  for 


COLUMNS 

FORMAT 

FORTRAN 

1-  9 

F9.4 

ST(K) 

10-18 

F9.4 

B(K) 

19-27 

F9.4 

H(K) 

28-36 

F9.4 

CA(K) 

37-45 

19 

NM(K) 

Explanation 

Any  identification  to  be  printed  at  the  top  of 
each  page  of  the  output 


Exp' ^nation 

L  =  length  between  perpendiculars  of  each  hull  (ft) 
Minimum  \ 

Maximum  1  nondimens ional  wave  frequency: 

Increment  of  ; 

Linear  ratio  for  converting  input  dimensions  on 
Card  Sets  3  and  4  to  the  size  ship  specified  in 
Columns  1-9  of  this  card 

n  =  number  of  stations  =  number  of  cards  in 
s 

Set  3 

Control  for  added  mass  (a)  and  damping 
coefficient  (b) 

If  IAMD=0,  a  and  b  will  be  calculated  by  this 
program  by  using  the  Lewis -form  method  for  the 
sections. 

If  IAMD=1,  values  of  a  and  b  are  input  in  Set  5. 
each  station) 

Explanation 

Sta^  =  station  number  ' 

b^  =  full  beam  at  the  waterline  \ 

h,  =  distance  from  keel  to  waterlinelomit  if 

U>0  c* 

C.  =  area  coefficient  =  r— — -r-  * 

\  bk  •  hk  >- 

ii 

m,  =  number  of  waterJines  for  which  offsets 
are  given  in  Card  Set  4 

If  values  of  b,  ,h,  ,  and  C.  are  given,  then 
m,  =0.  k  K  \ 


Note:  Cards  in  Set  3  must  be  in  order  of  ascending  station  numbers. 

Stations  0  (at  FP),  10  (at  amidships),  and  20  (at  AP)  must  be  in¬ 
cluded,  together  with  enough  additional  stations  to  define  the 
sectional  area  curve.  The  maximum  number  of  stations  i-'  30. 
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<»  in 


CARO  SET  4  (one  subset  for  each  station  with  m^  >  0) 

COLUMNS  FORMAT  FORTRAN  Explanation 

1-72*  8F9.4  Y(J,K)  y^  ^  ^=1,^)  =  half  beam  at  z.  ^ 

1-72*  8F9.4  Z(J,K)  z.  ^  (j^l.m^)  =  distance  above  the  baseline 

Note:  Values  of  z  must  be  in  ascending  order,  with  ^  at  the  keel  and 
z  ,  at  the  waterline.  * 

vk 

Subsets  must  be  in  order  of  ascending  station  numbers. 

If  m^  =  0,  there  will  be  no  cards  in  the  subset. 

If  >  0,  there  will  be  2,  4,  or  6  cards  in  the  subset. 


CARD  SET  5 

COLUMNS  FORMAT  FORTRAN 
1-72*  8F9.4  A33(N) 

1-72*  8F9.4  B33(N) 


Explanation 


a 

n 


(n=l,np)  =  added  mass  /  (pV) 
(n=l,np)  =  damping  coefficient 


Note:  np  =  number  of  wave  frequencies  =  ((OMAX-OMIN)/DOM)  +  1 

If  IAMD=0,  there  will  be  no  cards  in  this  set. 

If  IAMD=1,  the  values  of  a  and  b  can  be  obtained  from  the 

n  n 

columns  labeled  A33  and  B33,  respectively,  of  the  output  from 
Program  YF17  which  uses  either  the  Lewis-form  or  the  close-fit 
method  as  desired. 


CARD  SET  (  (one  card) 


COLUMNS 

FORMAT 

FORTRAN 

Explanation 

1-2 

12 

NSWH 

n^  =  number  of  significant  wave  heights  for 
irregular  sea  computations  <  4 

CARP  SET  7  (one 

:  cara) 

COLUMNS 

FORMAT 

FORTRAN 

Explanation 

1-36 

4F9.4 

H13(M) 

^1/3  =  significant  wave  height  (ft) 

•  m 

CARD  SET  8  (one  card) 

COLUMNS 

FORMAT 

FORTRAN 

Explanation 

1-2 

12 

NC 

n  =  number  of  conditions  =  number  of  cards  in 

c  Set  9 


* 

Continue  on  additional  cards  if  necessary. 
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CARD  SET  9  (one  card  for  each  condition) 


COLUMNS 

FORMAT 

FORTRAN 

Explanation 

1-  9 

F9.4 

BP 

L  = 

length  between  perpendiculars  (ft) 

10-18 

F9.4 

DK 

KD  = 

distance  of  deck  above  keel  (ft)  | 

19-27 

F9.4 

GK 

KG  = 

distance  of  CG  above  keel  (ft) 

,  at  LCB 

28-36 

F9.4 

BK 

KB  = 

distance  of  CB  above  keel  (ft) 

37-45 

F9.4 

YLP 

Bd  = 

transverse  distance  from  CL  of  catamaran 

u 

to  outer  edge  of  deck  (ft) 

46-54 

F9.4 

YL 

B  = 

transverse  distance  from  CL  of  catamaran  to 

c 

CL  of  one  hull  (ft) 

55-63 

F9.4 

RG 

V 

transverse  gyradius  (ft) 

Note:  The  irregular  sea  computations  are  done  for  the  ship  length  specified 
on  this  card.  This  L  can  differ  from  the  value  on  Card  1  which  is 
used  only  for  nondimensional  computations. 


APPENDIX  C 


SAMPLE  OUTPUT  FROM  PROGRAM  RLAC 
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Preceding  page  blank 


25 


9(«l  •  *e*  SPeCTRUM  IPICRSOM-HOSKOMITZ J 


VKATI CAL  ACCSLKMATIONS  ,Vt  A  CATAMARAN  IN  SCAM  SEAS  -  ASM  ASYMMETRICAL  MULL-rORM  -  MODEL  SCSI 


m  ooo<"«n«inn0»««0N»N4)nnMNMMOono<«^>«<Aooooooooooooooo 

«  ooooo**n4*«9M«M>>KN««oooooooooooooooooooooooooooooo 

<  oooooooo-*Mfi«nioooooooooooooooooooooeooooooooooo 

O  OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

. . . 

«  ©  o  o  o  o  o  o  o  o  o  o  o  ©  o  o  o  o  o  a  o  o  o  o  o  o  o  o  ©  o  o  o  9  o  o  o  oo  o  o  o  o  o  o  o  o  o  o  o  o  o 


-N  90K©w*VM#MO*Ko®N*«*nfii**i»»©«oi*>M©«rN««*#F»»'»N«i«MMM  —  MMM© 

b  oonA*Bn«B4«N<4n«nor«A«nNN^«^Moooooooooooooooooooo 

«cooooc*(Mi-OK»>*nw--.-oooooot)Oooooooooooooooooooooooo 


—  o 
►  o 
IL  O 


•  m 

X  o 


a  <f 


«  M 


I  M  OOOOOnNB^««>A«09tnnNMM**«*0000*A<4MMOOOOOOOOOOOOOOO 

«  •  oooooo-«n9>®co«Noooooe>ooooQ^9oooooooooooooooooo 

I  <  PPOOPPPPOmmmPOOOOOOOOOPOOOOOPOOOPPPPOOPPOPPPO 

I  o  POOPOOOOOOPOOOPOOOOOOOOOOPOOOOOOOOOOOOOOOOOOO 

•  I  «•  t  •«•••  i  •(  t  r  •••  t  ••  t  •••  #  ««•«•*•••  A  •••«•  t  ••«••••••>  S 

►  I  ffppppppppopppppppppppopppppppopoopppppopppppopopppp 

B.  I 

P  I  M  OOOO-WOl'INK^WWf'l-OOOOOOOOOOOOOOOOOOOOOOOQOOOOO 

•  I  •  OOPOOO«4P>PPKP)««<fOOOOOOOOOOOOOOPOPPOapOOOPOOOOOP 

« i |  opppopopoopoooooooo^ooppoooooppppepppooppoopp 

I O  OPPOPPPPOPPPPOOOOOPPOPPOOOPOPPPOPPPPOCOPOPPPP 

I  4  t  •  I  ••••••••••*••••  #  «••••••«•••  f  t  ••  t  t  . 

■  s  (tooppoopopooppopooooooooooopoocaoopooooaoooooooppoo 

I 

«  I 

Ml  OPPOONOpN»PO9t>'>Ofr«P*’llP)<VpM0O#MpP<V«Otf)  — 

Z  I  -  OaPOMOPK0m«®N00NNm0otnMK0MMP®K0#>««*)’’tft|<««MMMMMMMO 

i  b  eooppM^KowMMppKptfi^ninniAfMMMMpppooooooopppopopooP 


I  tflOPPOPOOOOPPpMMMMMPoaoPPOOPPOOPOOOPOOPOOPOOOOPPPOOP 
t 


0  BJP*.tf)P»M0C«®mM(B«4>«nM0>K««*lMpK0#A|0PK*>*)NO<aK*>mM00®«nM»00# 

X  *»  XPMM^in«0«KNPPO*  M<VP)««tf>«K0CPOMMNP)4rtf)00H0B9OM<*(V’*)«tf)0®fwpp 


NNNNNNNNM 


p  XMM«*MNNM")m«4>pN»opopnNM  ooooooooo  loooor^onooonoooonoo 

m  n  . . .  •  . . .  . . . . •  •  •  • 

XN  Z^mmmmmmmmmmmmmIBmmOOPOOPOOOPOPPOO^OOPCOOOPPOOOOOOOP 

-j  •  a 


.J®*>«tM®«M®«KO<»*OOM«®4rP««©K*n<VO»N'Otf)«#nNMMnO©»a0®KKK*’« 

P  »K#KPPP«|Immmm®POOOPOOPOOOOOOOOOOOOCOOOPPOOPOPOOOOPO 

M  p  fl  P|  M 

►  P  M 

B.  •  M 

w  O  poooopcpooooppoopoopoooooooppoooooppopooooooooooopo 

J  m  ,N««BO^##*ON4«aOM*«BON4«CON4«BON4«BOM#«JON««BON«4UiO 
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APPENDIX  D 

FORTRAN  LISTING  OF  PROGRAM  RLAC 


29  Preceding  page  blank 


CCMN0N/8L1/TITLE(12) 

Cfci-MUfl/eu2/NFR.UMLG<  50  >  •  ACMH(  50  I  .  D/MH( 50  ) 

DIMENSION  NM (  30  )  . 8 (  33  >  .  H (  30  >  . C A(  30  )  .  AR{  30  )  •  ST (  30  )  •  X (  30  )  . AM  (  30  )  . 

1  a3(30).SHB(30>.OS(3Q).Y'20.30).Z(23.30)*CM(50)*A33(50) .833(50). 

2  H13(4) <FRCLL(50 .4  i * FACC { 50 . 4 ) , RDEG ( A ) . AVGt  A ) * AMGt * ) * SM ( A ) • AD6G( 4 ) 
RA0*57.2958 

RHC* I *9905 

G*  32*1 74 

RHG=>RH04G 

P I *3 • 1415926 

PI2*2.4PI 

PI2G=PI2*G 

G81«.0081*G*G 

1  READ  (5.500)  (TITLE! J).J=1. 12) 

READ  (5.502)  BPL  .ON  1  N.CM  AX  j  DOM  .  CS7  <>  NOS  »  1  AMD 

READ  (5.504)  { ST (K ) . B( K ) .H ( K ) . C A ( K ) ,NM( K ) . X* 1 . NOS ) 

NUX* NM ( 1 ) 

DO  8  K=2.NCS 

IP  (ST(K).NE.IO.)  GO  TO  5 
MS— K 

S  IP  ( NUX .GE.NM(K) )  GO  TO  8 
NUX*NM(K) 

8  CONTINUE 

IP  (MS.NE.O)  GO  TO  10 
MR  ITE(6  <606 ) 

GO  TO  1 

10  SS*B PL/20 • 

DO  15  K=I.NOS 

IP  (NM(K).NE.O)  GO  TO  11 

B(K)=B(K)*CST 

H(K)«H(K)4CST 

AR(K)sCA(K)*e(K)4H(K) 

GO  TO  13 

11  N’zNM(K) 

READ  (5.506)  ( Y ( J. K ) . J= 1 , NZ ) 

READ  (5.506)  (Z( J.Kl.J-l.NZ) 

H(K)=Z(NZ.K)-Z(1 *K) 

B(K)*2.*V(NZ.K> 

AR(K)=2.4SIMPUN(Z(  1  .K)  ,Y(  ).  ,K).NZ> 

B ( K ) *B( K ) *CST 
H(K) =H(K)*CST 
AR(K)sAR(K) *CST 442 
CA(K)=AR(K)/(6(K)«H(K) ) 

13  X(K)=SS*ST(K> 

AM(K)=X(K) AAR ( K ) 

15  B3(K)i(CiS*E(K) )443 

IF  (IAMO.GT.0)  GO  TO  31 

16  00  30  K=1 .NOS 

IF  ( B ( K ) .LE.O . O.CR.H(K) »LE . 0 • )  GO  TO  30 
AC*CA(K) 

RAT=O.S*B(K)/H(K) 

TAR*1 ./RAT 

IF  (RAT.LE.l.)  GO  TO  23 
BL=0.29456*(2.-TAR) 

GO  TO  24 
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23  BL*0.294364C2.-RAT> 

24  UL«0*09812S*CRA' 4TAR410. ) 

IF  CCACKl.Gf.BL )  CO  TO  23 
CACK)*BL+0«300i 

GO  TG  26 

25  IFICACIO.LT. UL>  CO  TO  30 
CACK)*UL -0*0001 

26  HRITEC6.608:  ST C K > . AC *CA CK ) 

AR  CK ) *C  A  { K ) *8  C  K ) *H ( K ) 

AMK)*XCK)*ARCK  j 

30  CONTINUE 

31  CONTINUE 

HR ITE (6  *600 )  CTITLECJI.J31.12) 

HR I TEC 6. 602 ) 

00  33  K*1 . NOS 

35  HRITEC6.604)  ST C K) *BCK ) r HC K ) . ARC K ) .CA C K ) 
VOL*  5  I MPUN  CX.AR.NOS) 

BM*8CMS> 

HM*H  C  MS  > 

CB*VOL/CBPL*BH*HN) 

AH*S I MPUN C  X .B.NOS ) 

CH*A»/C 8PL4BH ) 

BOY*  S 1 MPUNC  X, AN, NOS) /VOL 
CBL*BOY/ePL 
OO  38  K*1.L0S 
38  SHECK)*CXCK)>80Y)4BCK) 

FLC*BCV*SIHPUNC  X.SHP.NQSI/AH 

CFL*FLC/6PL 

V0L=V0L*2. 

OLBS=RHG*VCL 

O CP* SI MPUN CX.B3. NOS) *0.6666667 

TM*VOL*RHQ 

BL*BPL/fiM 

8T*BM/HM 

OTCNS=DLOS/22  1. 

HR  I T  E  C  6*600 )  CTITLEC J). J*t*12) 

HR  ITE  C  6*610)  OPL.BM.HM.VOL.TM 
HRITEC6.6) 2!  DLBS.DTONS. AH.OIP.CBrCH 
HR  IT£(6*614)  BOY , CBL *FLC > CFL 
HA ITE  C  6  *C1 S ) 

HRITCC6.616) 

HR ITE (6 • 61 8 ) 

VCLNO=VOL/ C  2 . *3PL**3 ) 

00  40  K=1 • NOS 
BCK)*B( K)/EPL 
40  HCK)*M(X)/6PL 

NFR*  CCMAX-CM|N)/00M*I  ,2 
0MLGC1 )=OMIN 
B33C  l  >0. 

A33C  1)0* 

OQ  45  N-2  *  NFR 
OMLGCM  =CMLG(N-1 )*DOM 
A33C  N)=0. 

833C N)*0* 

45  CONTINUE 

IF  C’.AMD.GT.O)  GO  TO  56 
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Call  NILS(NOS.MS.ST.OS. JFK) 

IF  (JFK. EC. 0)  GQ  TO  1 
DO  50  K*1 .NOS 

CALL  ADM AB (8(KI «K(K) . CA ( K ) ) 

DO  50  N*1.NFR 

A33(N)*A33(N ) 4DS( K I */OKH( N ) 

50  B33(N)*B33(N)4DS(K>*0AKH(N) 

00  55  N*l.NFR 

A33(N)*A33(N)/( OMLM N ) 44 24 VOL NO ) 

55  B33(N)*B33(N)/(0MLG(N)4VCLND) 

GO  TO  56 

56  READ  (5.506)  ( A33( N ) .N* I .NFR ) 

READ  (5.506)  (BJ3(N). N* 1 . NFR ) 

56  HZL=0.02 
MS*H2L4PI 
READ  (5.510)  NSMH 
READ  (S.508)  { Mi 3( M ) . N*l . NSMH ) 

READ  (5*5X0)  NC 
DO  100  NCD*1 »NC 

READ  (5.508)  BP .OK  »GK. EK • YLP. VL  *  RG 
IF  (BP.EO.ePL)  GO  TO  60 
RL*BP/8PL 
BPL=BP 
RLZ'RLARL 
RL-VRL2*RL 
Rl «=RL3*RL 
V0L=V0L*RL3 
AM*AM4RLZ 
0 1 P=  OIPARLA 
GLOSSOLOGIES 
DTCNS*D  (0NS4RL3 
TM*TM4RL3 
BK=BM#RL 
HM*HM4RL 
OK=D  :^4RL 
8K=BK*RL 
60  VOLI =V0L/2. 

TKGL*TM4SQRT ( G/BPL ) 

YLB* YL/BM 
YL2=YL*4? 

SRLG-SQRT ( BPL/G ) 

GD*DK— GK 

01*2 .4(CIP4YL24AM) 

RG2=RG*RG 

RGL- RG/YL 

SEP* (2.*VL-BM)/BM 

TM I*RC2  4TM 

BG*GK  BK 

8MM*0I/V0L 

GM*BK4BMM-GK 

CRM* ( 01 /VOL-BG) 40LBS 

OFV=OIP—BG4VOL 1 

MR  JTE (6.600 )  (TITLE! J). J*I.I2) 

MR  ITE ( 6.620 )  BPL.8L.BT.SEP.VLP.0K.GK.BG.GM.RGL.0T0NS 
MR ITE (6 .623 )  ( HI 3( M ) . M=I , A ) 

00  80  N*I .NFR 
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ON(N)«OMLG(N)/SRLG 

0M2«0M(N)#*2 

OM4»OM2*OM2 

0MS«0M4*0M(N) 

ML*P I 2G/0M2 

«LL>  ML/BPL 

MH2»H2L*ML 

MH-0.S4MH2 

ANG«PI2*YL/ML 

PIL>PI2/«L 

FBAR*MH24RHG*SQRT(  (VL*AM4EXP(-PML*V0LI/AM)4SIN(ANC)  )*42  ♦ 

1  (PML4QFV4C0S( ANGl ) **2 I 
CMaTMI*TL2*A33(N)*TM 
CN*VL?.*833 ( N) 4TMGL 

PH  IB>lFBAR/SQRT  (  ( CRM-CM4CM2  ) 4424CN44240M2  ) 

PHlBsABSIPHlB) 

PH I* PM 10/ US 
RACR*(PHIB/MH)**2 
RACAsRA0R*CM4 
00  70  M*l*NSMH 

Sl!Ni><Sei/ON«tEXP(>33t96/IH13(IO**nOM4)  ) 

FRCLL (N,V) «RA0R*SM ( M ) 

70  F ACC  (N.M)*RA0A*SM(M) 

BO  MR ITE (6*624 )  OMLG( N> • MU. .PHI  . OM(N)*( SMI M  )  * FROLL ( N*M ) *PACC (N*M ) • 

1  M*1 .4) 

00  85  Mal.NSVH 

E2»S IMPUN( CM. FROLL ( 1 • M ) * NFR ) 

ROEG ( M ) «SCRT ( E2 ) *2 .0  ARAD 
E2*StMPUN(CM.  F ACC (  1  •  M  )  »  NFR ) 

SRE«SQRT(E2)*2*0 
A0EG(P)-SRE4RA0 
AVG(M|«YLP*SRE/G 
AHG(M)*GO  ASRE/G 
BS  COKTINUE 

MRITEC6  *600 )  (TITLE!  J)*.'*l*12) 

MR ITE  (  6*620  )  BPL.BL  .  BT  .  SEP.  YLP.DK  *GK*  BG.  GM.RGL  .OTONS 

MR ITE (6  *630 )  (H13(M) *M* 1 .4 ) . ( ROEGt P )• M*l *  4 > • ( ADEG( M ) *M*l *4 1  * VLP. 

1  ( AVG(M).M*1 *4) *G0. ( AH5(M).N*1.4> 

100  CGATINUE 
GO  TO  1 

500  FORMAT  (12A6) 

502  FORMAT  (SF9. 3*219) 

504  FORMAT  (4F9.4.I9) 

506  FORMAT  (8F9.4) 

508  FORMAT  (8F9.3) 

StO  FORMAT  ( 12) 

600  FORMAT  (59H1  VERTICAL  ACCELERATIONS  ON  A  CATAMARAN  IN  BEAM  SEA 

IS  -  •  1 2A6 ) 

602  FORMAT  (1H0/73H  STATION  BEAM  (FT)  DRAFT  (FT)  ARE 

I A  (FT2)  AREA  COEFF.  ) 

604  FORMAT  (5FI4.3) 

606  FORMAT  ( 1H0.22HSTAT 1CN  10*0  NOT  GIVEN  ) 

608  F-ORMAT  ( 1H0  •  I  OHS  TAT  ION  *  ,  F9.4  *6X*  30HAREA  COEFFICIENT  CHANGED  FROM 
1  • F10*  4  *2X  *2HT0*2X  *F10  *4 ) 

610  FORMAT  (49H0  RHO  *  MATER  DENSITY  *  1*9905  LB*-SCC2/FT4*  15X* 

1  48HG  *  ACCELERATION  OF  GRAVITY  a  32*174  FT/SEC2  //  6X* 
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2  49HL  *  LENGTH  BETWEEN  PERPENDICULARS  OP  EACH  HULL. 15X. IH«. 

3  F10.3.3H  FT  //  6X»  SOHB  »  FULL  BEAM  OF  EACH  HULL  AT  MIOSH.'PS  ( 

4STA. 10) .S4X.1H».F10.3.3H  FT  //  6X.37HH  *  ORAFT  ( WL  TO  KEEL)  AT 

SMISDH1PS.  27X,lHa,Fl0.3.3M  FT  //  6X.46HV0L  »  VOLUME  OF  WATER  01  SP 
6LACE0  BY  BOTH  HULLS . IBX. 1H>, F10. 0 .  4H  FT3  //  6X.  44HM  *  TOTAL  M 
TASS  OF  CATAMARAN  a  RHO  *  VOL . 20X • 1H».F 10*0* I IH  LB-SEC2/FT  ) 

612  FORMAT  (61H0  D  a  DISPLACEMENT  (CROSS  WEIGHT)  OF  CATAMARAN  a 

\  M  0  G.  9X.iHa,Fk0.0.  6H  LB  a,  F10.1.5H  TONS//  6X.3SHAW  «  WATER 
2PLANE  AREA  OF  EACH  HULL . 29X. 1 Ha, F 10 .0 • 4H  FT2  //  6X.  6SHIW  a  POME 
3NT  OF  INERTIA  OF  AW  W/RESPECT  TO  LONG.AXIS  OF  HULL  a.FlO.O. 

4  4H  FT4  //  6X •  5SHCB  ■  BLOCK  COEFFICIENT  OF  EACH  HULL  •  VCL/2/ 
5(L*B*H) .9X.IH*.F10.3//  6X.  SSHCW  *  WATERPLANE  COEFFICIENT  OF  EAC 
6H  HULL  *  AW/ ( L4B ) • 9X • 1H* .FI 0 • 3 ) 

614  FORMAT  (71H0  LC8  a  LONGITUDINAL  CENTER  OF  BUOYANCY  (DISTANCE 

1 AFT  OF  FT)  a.F10*0*6H  FT  *.F7.3.2H  L  //  71H  LCF  a  LONG I 

2TUDINAL  CENTER  OF  FLOTATION  (DISTANCE  AFT  OF  FP)  a.FtO.O. 

36H  FT  *  »F7 • 3  »  2H  L  ) 

615  FORMAT  (60H0  OK  a  VERTICAL  DISTANCE  FROM  KEEL  TO  DECK  AT  MID 

1 SHIPS  //  69H  8K  a  VERTICAL  OISTANCE  FROM  KEEL  TO  CENTER  OF 

2BUCYANCY  (  C  .8* )  ) 

616  FORMAT (6SH0  GK  a  VERTICAL  OISTANCE  FROM  KEEL  TO  CENTER  OF  GR 

1 AV IT Y  ( C.G . )  //6X • S4HEG  a  VERTICAL  DISTANCE  FROM  C.8.  TO  C.G.  « 

2  GK  -  BK  //  6X •  SOHMG  =  METACENTRIC  HEIGHT  »  BM  ♦  8K  -  GX 

3  //  6Xt  6SHL1  a  HORIZONTAL  OISTANCE  FROM  CL  OF  CATAMARAN  TO  C 
4L  OF  ONE  HULL  //  6X.  69HL2  *  HORIZONTAL  OISTANCE  FROM  CL  OF  CATA 
SMARAN  TO  OUTER  EOGE  OF  OECK  //  6X,  38HSEP  »  MULL  SEPARATION  a  (2 
6  •  LI)  -  B  //  6X.  26HRG  a  TRANSVERSE  GYRAOIUS  //  6X.  30HH13  a  SI 
7GN IF ICANT  WAVE  HEIGHT  //  6X.  18HWL  a  HAVE  LENGTH* 10X, 19HWH  a  «AV 
8E  AMPLITUDE. 10X.46HWS  a  HAVE  SLOPE  -  2  *  3.14  *  WH  /  WL  »  3.14/50) 

618  FORMAT  (71M0  W  a  HAVE  FREQUENCY  a  FREQUENCY  OF  ENCOUNTER 

1 ( BEAM  SEAS  ONLY)  //  6X.30HPH1  a  AMPLITUDE  OF  ROLL  ANGLE  //  6X. 

2  39HS(H)  a  SEA  SPECTRUM  ( P IER50N— MOSKOH I TZ //  6X, 

3  26HRA0R4S  a  (PHI/MH)SQ  *  S(W)  •  1SX. 

4  30HRA0A4S  a  ( PH I 4HAH/HH ) SQ  •  S(W>  ) 

620  FORMAT  (128H0  LIFT  I  L/B  B/H  SEP/B  L2 (FT ) 

1  OK (FT )  GK(FT)  BG(FT)  GM(FT)  RG/Ll 

2  O(TONS)  /  3F 10.3* 7F 12 .3 • F 1 4 • 1  ) 

622  FORMAT  ( 1 HO ) 

623  FORMAT  (1H0.30X.4(7X,5FH13  a.FS.l.AH  FT..4X)  /  31 X. 4( 2X . 23( 1H- ) )  / 
1131H  W ( L/G )  WL/L  PHI/HS  HI  1/SEC)  SIW)  RA0R4S  RAOA4S  S( 
2H‘  RA0R4S  RA0A4S  SIH)  RA0R4S  RA0A4S  SIW)  RA0R4S  RAO 
3A4S  ) 

624  FORMAT ( F6.2 »F7» 2.F8. 2.F9.3.1X«4(F9.3. 2F8.4 ) ) 

630  FORMAT  t 1H0 *5X . 28HS I GN IF IC ANT  WAVE  HEIGHT  (F T ) .21 X. 4F10.2  / 

1  1 HO .SX * 32HSIGN IF ICANT  ROLL  ANGLE  ( OEGREES ) • 1 7X. 4F10.2  / 

2  l HO .SX.40HS1GNIF ICANT  ROLL  ACCELERATION  (0EG/SEC2) .9X.4F10.2/ 

3  1H0.5X.26HSIGN. VERTICAL  ACC.  /  G  (  .F6.1.12H  FT  FROM  CL).  SX 

4  .4F10.3  /  1H0.SX.26HSIGN.  HORIZ*  ACC.  /  G  (  .  F6.1. 

5  17H  FT  ABOVE  VCG)  .  4P10.3) 

ENO 
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